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We used a torsion pendulum and rotating attractor with 20-pole electron-spin distributions to 
probe dipole-dipole interactions mediated by exotic pseudo-Goldstone bosons with m\, < 500 peV 
and coupling strengths up to 14 orders of magnitude weaker than electromagnetism. This corre¬ 
sponds to symmetry-breaking scales F < 70 TeV, the highest reached in any laboratory experiment. 

We used an attractor with a 20-pole unpolarized mass distribution to improve laboratory bounds 
on CP-violating monopole-dipole interactions with 1.5 peV< mb < 400 peV by a factor of up to 
1000 . 
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Spontaneously-broken global symmetries play an im¬ 
portant role in particle physics [ 2 ]. When the underlying 
symmetry is exact, the process always produces mass¬ 
less pseudoscalar Goldstone bosons whose coupling to a 
fermion with mass to/ is g p = rrif/F, where F is the 
energy scale of the spontaneously broken symmetry. If 
the symmetry is not exact but explicitly broken as well, 
as in the chiral symmetry of QCD, the fermionic cou¬ 
plings are unchanged, but the resulting pseudo-Goldstone 
bosons, such as the QCD pions, acquire a small mass 
TOb = A 2 /F where A is the explicit symmetry-breaking 
scale of the effective Lagrangian. Searches for the ultra- 
weak, long-range interactions mediated by exotic pseudo- 
Goldstone bosons, therefore, provide very sensitive and 
general probes for new hidden symmetries broken at ex¬ 
tremely high energies. 

The tree-level potentials from pseudoscalar boson ex¬ 
change are purely spin-dependent. The classic pseu¬ 
doscalar potential is the dipole-dipole interaction 
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where A = h,/(in\,c). Axion-like bosons with an addi¬ 
tional scalar coupling, gs, can also generate a monopole- 
dipole interaction [2J 
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Because these potentials average to zero for unpolarized 
bodies, traditional searches for new macroscopic forces 
are essentially insensitive to such bosons. Motivated by 
theoretical conjectures that propose additional pseudo- 
Goldstone bosons such as axions, familons, majorons, 
closed-string axions and accidental pseudo-Goldstone 
bosons (see !3| for a recent review), we developed a 
generic “pseudo-Goldstone detector” with high sensitiv¬ 
ity to both Vdd and E m d interactions. We combined the 


strategies of previous Eot-Wash torsion-balance probes 
of electron-spin-dependent forces [3, [H (closed magnetic 
circuits containing high and low spin-density materi¬ 
als) and short-distance gravity 0 (planar geome¬ 
try, high-multipolarity signals and continuously rotating 
attractors) to produce the torsion-pendulum/rotating- 
attractor instrument shown in Fig. |T| The small scale 
of our device allowed us to probe Vdd interactions 
with TOb < 500 peV. Previous studies with polarized 
electrons [Sj and neutrons @, [§] were restricted to TOb ^ 
5 peV. 

The key element of our instrument was a spin-ring 





FIG. 1: Left: the 20-pole spin pendulum and spin-attractor; 
p-metal cans on the pendulum and attractor are cut away to 
show the Alnico (green) and SmCos (blue) segments and one 
of the 4 pairs of calibration cylinders (red). The mirror cube 
allowed us to monitor the pendulum twist angle. The mag¬ 
netic shield surrounding the entire pendulum is not shown. 
Lower and upper right: top views of the spin and mass at¬ 
tractors, respectively. Arrows indicate net spin density and 
direction. 
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containing 20 equally-magnetized segments of alternat¬ 
ing high and low spin-density materials. This formed 
a spin 20-pole with a negligible external magnetic field; 
the sensitivity to Vdd or Vmd interactions arises entirely 
from the spin-density contrast in the rings. One spin-ring 
was the active element of our detector, a torsion pendu¬ 
lum placed just above an attractor rotating at frequency 
w. Our dipole-dipole search used an attractor consist¬ 
ing of a second 20-pole spin-ring so that Vdd interactions 
would produce a 10u; torque on the pendulum. We used 
an unpolarized copper attractor that formed a mass 20- 
pole to measure the gravitational background in our Vdd 
study and to probe V m d interactions. The lOw torques 
from Vmd and from gravity were distinguishable because 
their turntable phases differ by 9 degrees. The pendulum 
and both attractors each contained four 4.8 mm-diameter 
cylinders (tungsten and vacuum for the spin and mass 
attractors, respectively) that provided a continuous 4w 
gravitational calibration signals. 

Alternating wedges of SmCos and Alnico provided the 
spin contrast of the rings. SmCos has a substantial 
orbital contribution to its magnetic field [3] while Al- 
nico’s magnetism comes almost entirely from polarized 
electrons. SmCos fully magnetized to 9.8 kG contains 
~ 4.5 • 10 22 spins/cnr 3 while Alnico magnetized to the 
same degree has ~ 8 • 10 22 spins/cm 3 [4J|. We cut the 
SmCos wedges from commercially magnetized material 
while the Alnico was magnetized simply by assembling 
the ring. We tuned the precise magnetization of each Al¬ 
nico wedge in situ by applying a localized external field 
until the peak-to-peak leakage field 3 mm from the ring 
was reduced from ~100 G to ~8 G. We then enclosed the 
ring assemblies in nested 2-layer /x-metal cans with a total 
thickness of 0.53 mm, reducing the peak-to-peak residual 
field to ~ 10/xG. A 0.99 mm thick shielding screen con¬ 
sisting of 21 layers of alternating /x-metal and aluminum 
foils separated the attractor assembly from the pendu¬ 
lum. A 1.27 mm thick, cylindrical, /x-metal “house” sur¬ 
rounded the pendulum except for a hole for the suspen¬ 
sion fiber and another that provided optical access to 
the pendulum; a 0.76mm-thick /x-metal tube surrounded 
the attractor turntable. The peak-to-peak magnetic field 
change at the pendulum location with the full shielding 
in place was below our resolution of 2 /xG. 

Could the magnetic shielding also shield the Vdd and 
Vmd interactions? We discuss this important point 
elsewhere[lO] and show that such a shielding effect is 
completely negligible in our case. 

The mass-density difference between SmCos (8.31 
g/cm 3 ) and Alnico (7.31 g/cm 3 ) would produce a sig¬ 
nificant gravitational 10w torque. We placed 76 /xm thick 
W shims above and below each Alnico segment and Ti 
shims above and below each SmCos segment to minimize 
10w gravitational torques. 

The active elements of the pendulum and attractor, 
along with the magnetic shielding, were installed in a 


rotating-attractor torsion balance normally used to study 
short-distance gravity. Details of that torsion balance 
and the general methods of data analysis are given in 
Refs.jl, 0. The only other significant change to the ap¬ 
paratus described in Ref. 0] was an improved attractor 
drive that locked the output of a 2 21 pulses/rev angle- 
encoder to a crystal-controlled oscillator[ll|. 

We centered and leveled the attractor ring to ~ 15 /xm 
of the turntable rotation axis with optical and mechani¬ 
cal techniques. We leveled the pendulum to 130 /xm of its 
rotation axis using capacitive techniques, and centered it 
to about ±20/xm of the attractor rotation axis by max¬ 
imizing the 4 u> signal from the mass attractor. These 
misalignments were negligible relative to our typical sep¬ 
arations s > 2 mm. We inferred the vertical separation, 
s, between the bottom of the pendulum ring magnets 
and the top of the attractor magnets (or copper) using a 
^-micrometer on the vertical translation stage that sup¬ 
ported the suspension fiber. We measured the attractor- 
pendulum capacitance as a function of z and fitted these 
data to a finite-element electrostatic model to map the 
^-micrometer readings into the pendulum-screen separa¬ 
tion. Mechanical and optical measurements provided the 
additional information needed to determine s. 

The torque on the pendulum was inferred from a har¬ 
monic analysis of its twist angle, 6 , as a function of the 
attractor angle <fi = Lot. The data analysis procedure was 
similar to that used in Ref. 0 ]. The 0(0) time series 
was processed by a 4-point digital filter that suppressed 
free torsional oscillations as well as the DC response and 
linear drift. Data runs were divided into cuts containing 
exactly one attractor revolution and each cut was fitted 
with a quadratic drift term plus the first 14 harmonics 
of the turntable angle. The harmonic amplitudes were 
corrected for pendulum inertia, electronic time constants 
and the response of the digital filter, and then converted 
into torques using the effective value of the fiber’s tor¬ 
sional constant n = I(2n/To) 2 = 3.1 aNm/nrad, where 
the To is pendulum’s free oscillation period and its mo¬ 
ment of inertia, I = 134 gem 2 , was computed from 
a detailed numerical model. To was determined from 
“sweep runs” taken after each science run; the attrac¬ 
tor turntable was stopped and the pendulum was given 
a ~ 10/xrad kick. The resulting oscillations were ana¬ 
lyzed to obtain a precise period and, crucially, to map 
out small nonlinearities in the autocollimator’s analog 
position-sensitive detector. The measured lOw and 4 lu 
torques from a run were found by weighting equally all 
(typically 48) cuts in that run with the statistical uncer¬ 
tainties determined by the scatter of the results. These 
data were compared to those expected from Vdd and V„id 
interactions and from gravity. The expected torques, as¬ 
suming that the pendulum was aligned with the attrac¬ 
tor, were computed using the Fourier-Bessel expansion 
which converges rapidly for our application. 

The attractor rotation periods, T at t = 7 Tq or 6To, were 
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FIG. 2: Sample power spectral densities of the twist signals 
from the spin and mass attractors at the closest attained 
separations. The dashed lines show the thermal noise, solid 
lines include the effect of an additional l// 2 component. The 
8ui peak in the mass attractor data is the 1st harmonic of the 
calibration signal. It is much smaller in the spin attractor 
data because of its larger value of s. 


selected to place our 4uj and 10w signals in low-noise re¬ 
gions (see Fig. [2|. The noise was dominated by ther¬ 
mal fluctuations from internal losses in the suspension 
fiber which gave the torsion oscillator a quality factor of 
Q « 1500. 

Spin-attractor data were taken at s = 4.12, 5.13, and 
8.15 mm (uncertainties are ±0.015 mm). Because of the 
1/r 3 fall off of the potential, our Vdd sensitivity comes 
entirely from the s = 4.12 mm data. The results from 
165 hours of s = 4.12 mm data are shown in Tabled] 
We expected the largest systematic effects with the spin 
attractor to be residual gravitational and magnetic cou¬ 
plings betweens the pendulum and the attractor. Mass 
attractor data supplemented by calculations showed that 
our shims reduced the gravitational component of Aio^ 
by two orders of magnitude to ~laNm. Measurements 
showed that the magnetic leakage field was fairly con¬ 
stant across all higher (> 5 u>) harmonics. We observed 
little evidence for such couplings (see Fig. [3]), which 
would have produced torques at all these frequencies as 


TABLE I: Observed 4ta and 10a; torques. Amplitudes A, are 
in units of aN m, phases (f> are in degrees, and separations s 
are in mm. The lcr uncertainties do not include systematic 
effects. If Fmd = 0, we expect A <f> = (pioui — <j>4ui = —9.0°. 
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FIG. 3: Comparison of the spin-attractor 10u; science signal 
with nearby background signals. The shaded horizontal band 
indicates the mean and standard deviation, er, of the back¬ 
ground signals. The horizontal line shows the mean ampli¬ 
tude, fry / 7r/2, expected for random signals whose quadrature 
components have zero mean and spread a. 

well as at 10w. As a result, no corrections for magnetic 
backgrounds were necessary. Other systematic concerns 
such as thermal and electrostatic effects were found to 
be negligible. Our final value is Aio u = (1.3 ± 2.2) aN nr 
for a 2cr upper limit of 5.1 aN m. The corresponding con¬ 
straints on a new dipole-dipole interaction, and the as¬ 
sociated bounds on the symmetry-breaking scale F, are 
shown in Fig 2] These are the most sensitive labora¬ 
tory constraints on (g®) 2 /?ic for m b < 500^eV/c 2 (at the 
5.5 x 10~ 17 level for mb < 30 /xeV). To our knowledge, the 
only other laboratory constraints on pseudoscalars in this 
mass range are Ramsey’s 1979 limit, ( g?) 2 /Tic < 3 x 10~ 4 
level [13, on anomalous spin-spin interactions between 
protons. Our results indicate that F > 70 TeV. 

Because magnetic backgrounds in our Vmd study us¬ 
ing the mass attractor were small, we could use a single 
0.25 mm thick /z-metal screen. This allowed us to take 
mass-attractor data at s = 1.98 mm as well as at 2.03, 
3.00, 4.04 and 7.99 mm (uncertainties are ±.015 mm). 
Our V m d constraints come entirely from the s = 1.98 mm 
data. The other data allowed us to check for systematics 
and validate our gravitational calculations. The observed 
4 ui and lOw torques from 38 hours of s = 1.98 mm data 
are shown in Table U] The gravitational contribution to 
the 10w torque vanishes when the copper arms of the at¬ 
tractor are directly below either the SmCos or the Alnico. 
This occurs (see Fig. [TJ when the attractor is rotated 9° 
away from the angle at which calibration cylinders are 
aligned. Conversely, the Vmd torque is maximal at those 
orientations. This allowed us to separate the effects of 
gravity from a V m d interaction. The V m d component of 
the 10 lj torque is 

Amd = Ai Ow |sinlO(A0±5</>)| , (3) 

where A <f> = (f> iq w — (j>/uo and 5<j) is nominally 9°. The 
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FIG. 4: Bottom: exotic dipole-dipole limits from this work 
and Ref. 5|. Arrows indicate the infinite-range constraints 
from Refs, [l^. m. Electron g —2 constraints are at the 10 10 
level [ 13 ], Top: limits on the symmetry-breaking scale from 
this work and Refs. 0,0. The shaded areas are excluded 
at 2<t. 


absolute value occurs because of the 4-fold ambiguity in 
the attractor angle inferred from the 4w signal. Our A m( j 
bound is dominated by the systematic uncertainty in 5<p. 
Alignment microscope measurements showed that phase 
of the magnet ring relative to the calibration cylinders 
was only fixed to ±0.17°. An estimated 50/im accuracy 



FIG. 5: Monopole-dipole constraints from this work and 
refs. [l3 - ^H The shaded region is excluded at 2o\ The mb = 0 
limit from Ref. 0]is 2x 10 -36 . (We doubled the la limits given 
in refs.0, HJ].) 


in positioning the gravitational shims, revealed by the be¬ 
havior of <^>ioco in our centering data, contributed an addi¬ 
tional error of ±0.29° and increased the uncertainty in Scj) 
to ±0.34°. This gives a ler result A md = (18 ± 12) aNm 
with a 2cr limit, A m( j < 38aNm. Our Kg^g^l/ftc con¬ 
straint, shown in Fig. [5] improves upon previous work by 
up to a factor of 1000 for 1.5 < mb < 400 /reV/c 2 . The 
most sensitive limit on (g^g^)/Tic is also at the 10 -28 
level [23 - 

Stellar cooling rates (23| constrain Vdd interactions of 
simple pseudoscalar particles at a level well below our 
bound, and the astrophysics bound on g®, combined with 
bounds on g^ from gravitational experiments, set very 
tight limits on Vmd interactions between electrons and 
nucleons[2J|. However, a chameleon mechanism could 
invalidate these astrophysical bounds while having a neg¬ 
ligible effect in cooler, less dense lab environments 25J. In 
this case Vdd and Vmd can only be constrained by labora¬ 
tory experiments such as this work which reveals that any 
hidden symmetry involving electrons must be broken at 
an energy scale F > 70 TeV and, if it is explicitly broken 
as well, that scale A must be > 0.1 MeV. These set the 
highest laboratory bounds on the minimum energy scale 
of new hidden symmetries involving leptons. Extensions 
of general relativity that include torsion as well as curva¬ 
ture predict infinite-range dipole-dipole interactions [26| 
and are also constrained by this work. 
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